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ABSTRACT 

We present high-resolution Keck optical spectra of the very young substellar eclipsing binary 
2MASS J05352184— 0546085, obtained during eclipse of the lower-mass (secondary) brown dwarf. 
The observations yield the spectrum of the higher-mass (primary) brown dwarf alone, with neg- 
ligible (~1.6%) contamination by the secondary. We perform a simultaneous fine-analysis of the 
TiO-e band and the red lobe of the K I doublet, using state-of-the-art PHOENIX DUSTY and 
COND synthetic spectra. Comparing the effective temperature and surface gravity derived from 
these fits to the empirically determined surface gravity of the primary (logg = 3.5) then allows 
us to test the model spectra as well as probe the prevailing photospheric conditions. We find 
that: (1) fits to TiO-e alone imply T cff = 2500±50K; (2) at this T off , fits to K I imply log g = 
3.0, 0.5 dex lower than the true value; and (3) at the true log g, K I fits yield T e s = 2650±50K, 
^150K higher than from TiO-e alone. On the one hand, these are the trends expected in the 
presence of cool spots covering a large fraction of the primary's surface (as theorized previously 
to explain the observed T e g- reversal between the primary and secondary). Specifically, our re- 
sults can be reproduced by an unspotted stellar photosphere with T e g — 2700K and (empirical) 
log g = 3.5, coupled with axisymmetric cool spots that are 15% cooler (2300K), have an effective 
logg = 3.0 (0.5 dex lower than photospheric), and cover 70% of the surface. On the other hand, 
the trends in our analysis can also be reproduced by model opacity errors: there are lacks in 
the synthetic TiO-e opacities, at least for higher-gravity field dwarfs. Stringently discriminating 
between the two possibilities requires combining the present results with an equivalent analysis 
of the secondary (predicted to be relatively unspotted compared to the primary). 

Subject headings: binaries: eclipsing - stars: low-mass, brown dwarfs - stars: pre-main sequence 
- circumstellar matter - stars: fundamental parameters - techniques: spectroscopic 



1. Introduction 

2MASS J05352184-0546085 (henceforth 2 M0535), a very youn g system located in the Orion Nebula 
Cluster (ONC), has recently been identified by Stassun et al. ( 20061 hereafter SMV06) as the first known 



substellar eclipsing binary (EB). EBs allow exquisitely precise direct determinations of the component masses 
and radii, and thus the surface gravities (log g), as well as the ratio of their luminosities (or equivalently, ratio 
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of their effective temperatures, T e g). As such, 2M0535 permits the first stringent tests of both the theoretical 
evolutionary models and the synthetic spectra that are widely employed to characterize the vast majority of 
brown dwarfs (for which direct measurements of mass, radius, and surface gravity are not possible). 



The parameters of the system found b y SMV06 w ere refined with more data bv lStassun et al.l (|2007l here- 
after SMV07) and still further bv lGomez Maaueo Chew et al.l (|2009l . hereafter G09). SM V07 found a spec- 
tral t ype of M6.5±0.5 for the primary (higher-mass) component, suggesting T e g w 2700K (jGolimowski et al 



20041 ). The latest analysis by G09 confirms the initial findings that: (1) Both components of 2M0535 are 
moderate mass brown dwarfs (Mi = 0.0572±0.0033 M Q , M 2 = 0.0366±0.0022 M ); (2) their radii {R x = 
0.690±0.011 R©, i?2 = 0.540±0.009 Rq) are consistent with the theoretical prediction that young brown 
dwarfs of a given mass should be much larger than their field counterparts^]; and (3) the T e g ratio of the 
components (T e s,2/T e s,i — 1.050±0.004) shows an unexpected reversal, with the primary being cooler than 
the lower-mass secondary. 

The reversal i n tem peratures is not predicted by any set of theoretical evolutionary tracks. To explain 



it. IChabrier et al.l (|2007l hereafter CGB07) proposed that strong magnetic fields on the primary suppress its 
interior convection and also produce cool surfac e spots; neither effec t is included in the standard evolutionary 
models, and both would act to depress its T e g. iReiners et al.l (|2007l hereafter R07) subsequently found that, 
compared to the secondary, the primary is a relatively fast rotato r with strong chromospheric Ha emission, 
which supports the presence of strong magnetic fields in the latter. G09l then showed that the observed small- 
amplitude residual (non-eclipse) variations in the 2M0535 lightcurve, modulated at the rotational periods of 
the primary and secondary, can be well reproduced by cool spots asymmetrically covering a small fraction 
(5, 1 0%) of both components' surfaces. While such small spots cannot explain the temperature reversal, 
Goal's analysis does indicate that spots are at least present. Moreover, they cannot rule out the very large 
(~ 50% areal coverage) spots on the primary required to explain the temperature reversal, as long as these 
are arranged symmetrically about the rotation axis (e.g., polar spots, equatorial bands, or 'leopard spots'). 



Additionally, through analysis of the optical to mid-IR spectral energy distribution (SED). lMohantv et al 



(|2009L hereafte r MS M09) showed that: ongoing accretion is highly unlikely in the 2M0535 system, lending 
credence to the lR07l conclusion that the Ha emission in the primary is chromospheric; and the system SED is 
consistent with effective temperatures of [T s,i, Tgg^] ~ [2700, 2900]K, in agreement with t he mid-M spec- 
tral types of the components rather than with the much lower ~ [2300, 2450]K proposed by |CGB07J within 
their the ory of T e g reversal via magnetic field effects. Combining their results with the others cited here, 
MSM09I concluded that while magnetically-induced spot /convection effects probably do play an important 
role in determining the T e g of the 2M0535 prim ary, as advocated by ICGB07I (and as indeed seems to be the 
case for active field dwarfs: Morales et al.ll2008l ). the theory is as yet insufficiently developed quantitatively, 
and small age variations between the components may play a significant role as well (as in fact appe ars to 
be the case in another young EB, albeit with stellar-mass components, Par 1802; IStassun et al.ll2008l ). 



Finally, subsequent to MSM09's work, MacDonald & Mullan (2009) have proposed a theory wherein 
magnetic fields inhibit the onset of convection (though do not suppress it entirely) throughout the 2M0535 
primary, instead of just in the upper-most super-adiabatic layers as in the theory of CGB07. The theory 
appears to reproduce the observations of 2M0535AB without invoking non-coevality or very large surface 
spots on the primary; this competing scenario must be evaluated as well. 



1 The 2M0535 component rad ii are slightly — ~10% — underpredi cted by some models fsee lSMV07ri ; the discrepancy becomes 
slightly stronger using the latest G09 radii cited here compared to SMV07's numbers. 
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In this paper, we present Keck HIRES observations of 2M0535 obtained during eclipse so as to isolate 
the spectrum of the higher-mass, lower- T c ff primary brown dwarf. Comparing the observed spectrum with 
state-of-the-art brown-dwarf atmosphere models, we test t he ab ility of these models to correctly reproduce 
the accurately known surface gravity (log 5 = 3.52±0.03; [gO^), and in the process we directly probe the 
prevailing photospheric conditions of the primary brown dwarf in the 2M0535 system. 



2. Observations and Data Reduction 

We observed 2M0535 on the night of UT 2007 Oct 23 with the High Resolution Echelle Spectrometer 
(HIRES) on Keck-fl We observed in the spectrograph's "red" (HIRESr) configuration with an echelle angle 
of —0.403 deg and a cross-disperser angle of 1.703 deg. In this configuration, the two features of primary 
interest in this paper, TiO AA8435-8455 and K I A7700, fall on the "green" chip, in echelle orders 42 and 
46, respectively. We used the OG530 order-blocking filter and the l'.'15x7'.'0 slit, and binned the chip during 
readout by 2 pixels in the dispersion direction. The resulting resolving power is R « 34 000, with a 3.7-pixel 
(~ 8.8 km s _1 ) FWHM resolution element. 

We obtained three consecutive integrations of 2M0535, each of 2400 s. ThAr arc lamp calibration 
exposures were obtained before and after the 2M0535 exposures, and sequences of bias and dome fiat- 
field exposures were obtained at the end of the night. The 2M0535 exposures were processed along with 
these calibrations using standard IRAF^ tasks and the MAKEE reduction package written for HIRES by 
T. Barlow. The latter includes optimal extraction of the orders as well as subtraction of the adjacent sky 
background. The three exposures of 2M0535 were processed separately and then median combined with 
cosmic-ray rejection into a single final spectrum. The signal-to-noise (S/N) of the final spectrum is ~ 15 per 
resolution element. 

Importantly, we intentionally chose the observations to coincide exactly with the secondary eclipse, i.e. 
when the lower-mass, smaller, higher-T c g secondary component was behind the primary as seen from Earth. 
The first exposure started at UT 12:20 h, and the third exposure ended at UT 14:22 h, corresponding to 
orbital phases of 0.0 709 and 0.0794, respectively, during which time the secondary is almost completely 



blocked (cf. Fig. 3 in[SMV07|). Integrated over the entire 2-hr observation, the total light contribution from 
the secondary was rs1.6%. The light contribution from the secondary was calculated using the accurately 
determined radius ratio, temperature ratio, and orbital parameters, including the orbital inclination, from 
the light curve modeling performed in G09. Thus the resulting spectrum is effectively that of the primary 
alone. 



3. Synthetic Spectra 

We use the latest version of synthetic spectra for plane-parallel atmospheres g enerated using th e 



PHOENIX code, designated AMES-Cond (version 2.4) and AMES-Dusty (version 2.4) (|Allard et al.ll200ll) . 
These synthetic spectra have become broadly used in the literature on low-mass stars and brown dwarfs 
especially at young ages. In addition, these model spectra are incorporated into the commonly used stellar 



2 Time allocation through NOAO via the NSF's Telescope System Instrumentation Program (TSIP). 

3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
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evolution models of iBaraffe et al.l (J1998) as well as in the CGB07 models discussed in Sec. Q] Thus one of 
our aims in selecting the PHOENIX synthetic spectra here is to assess these commonly used models in the 
context of brown dwarf evolution. 



The PHOENIX code jHauschildt fc Baron|[l999h is a general purpose stellar atmosphere model tool 
that makes use of very complex atomic models and line blanketing by hundreds of millions of atomic and 
molecular lines. The PHOENIX mo dels used h ere incorporate the most re cent AMES line lists for both 



TiO dLanghofj|l997i ISchwenke Ill998l) and H 2 (jPartridge fc Schwenkel[l997l ). A good treatment of H 2 is 
essential for analyzing optical spectra, even though H 2 opacity dominates only in the infrared (TiO opacity 
is more important in the optical). This is because the overall H 2 opacity is larger, and its lines occur closer 
to the peak of the SED than those of TiO, at the low T e g in M spectral types. Consequently, changes in the 
H 2 opacity have a substantial effect on the atmospheric temperature structure and thus on the emergent 
spectrum even in the optical. A total of about 500 million molecular lines are currently included in th e 
models; of these, -207 million are lines of H 2 0, and -172 million are of TiO (jAllard et aljboool boOll) . 
Here we use solar-metallicity models ([M/H]=0.0). While the metallicity of 2M0535 is not explicitly known, 
a large deviation from solar is not expected for a young object in a nearby star-forming region. We discuss 
potential metallicity effects on our results in more detail in §6.2. 

Dust formation is another potentially important effect in the low-temperature atmospheres of M-type 
objects; grains affect the atmospheric structure as well as the emergent flux. Bot h models we exam ine treat 
grain formation self-consistently, through chemical equilibrium calculations (see lAllard et al.ll2001l ). Under 
physical conditions where the chemical equations imply no grain formation, the COND and dusty spectra 
are identical; in the models, this occurs for T e gf > 2500K. For the latter temperatures, therefore, either set 
of synthetic spectra may be used. The difference between the two is in their treatment of dust settling, once 
grains have formed. The two models represent the two limits of settling: dusty models treat the case where 
grains form and remain suspended in the photosphere, while the COND ones are applicable when dust has 
formed but subsequently settled ( "condensed" ) out of the photosphere en tirely. Observatio ns of field dwarfs 



indicate that dust settling becomes important only in the L types (e.g.. ISchweitzerll200lh . For the mid-M 
spectral type of 2M0535, dust formation may occur but the grains are likely to remain in the photosphere 
I Jones fc Tsuiilll997i[MM ). We thus use dusty models for T eS < 2500K, and COND models for T cff > 2500K 
(where no dust formsjS 



4. Methodology 



We wish to determine the T e g and surface gravity (log g) of the higher- mass component of 2M 0535 (the 
"primary," hereafter 2M0535A) from comparisons to synthetic spectra. As iMohantv et al.l (|2004l hereafter 
M04) have shown, two ideal regions for this analysis are the TiO-e bandheads at AAA8435, 8445, 8455, and 
the red-lobe of the K I doublet at A7700 (the blue lobe falls in the gap between echelle orders in the HIRES 
setting used). In particular, the TiO bandheads are very sensitive to T e s, but negligibly so to gravity, while 
the K I absorption is sensitive to both; using the two regions in tandem therefore enables one to disentangle 
and individually determine these two parameters. 



Comparing the data to models requires some modifications to both. These are discussed in detail in 



4 As noted, DUSTY and COND are identical at these T e fj; at high-resolution, only COND models are available for T a g > 2500K, 
so we use them for T e ff > 2500K after verifying that they are indeed nearly identical to DUSTY at the overlap T e ff of 2500K. 
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M04c the salient points are as follows. The models are rotationally broadened (using Gray's methodology 



(Gray 1992), with a standard limb-darkening parameter of 0.6) by 10 km s 1 to match the observed v sini 



of 2M0535 A (jROTT ) , and further broadened by convolution with a Gaussian profile to match the instrumental 
broadening (finite resolution) of the data. Since our data are not flux-calibrated, comparison to the models 
also requires some form of scaling. This is accomplished by normalizing both the data and models by their 
average flux over a narrow region of pseudo-continuun[f| just outside the absorption lines of interest: over 
AA[8402. 5-8411.5] for the the TiO-e region, and over AA[7707.5-7709.5] for the K I region (wavelengths in 
the laboratory rest- frame). Recall that the data are also flat-fielded, which removes the blaze-function but 
preserves the innate shape of the stellar spectrum. Our normalization procedure then ensures that the data 
and models are only 'anchored' over a narrow wavelength range, but otherwise unconstrained, so the models 
need to match not only the absorption bands, but also the shape and slope of the continuum, to ensure a 
good fit. This provides an additional check on the veracity of the preferred fits. 

For comparisons to the stellar spectrum, we first use the DUSTY and COND models depending on the T e g 
being tested, as described in $3](cOND for > 2500K); the results are described in ij5] We then model the stellar 
spectrum as a combination of a naked photosphere and cool spots, as discussed in ij6.3l For this analysis, 
the photosphere and spot are represented by different spectra, depending on the adopted temperature of 
each (e.g., for a 2700K photosphere and a 2300K cool spot, we use COND for the former and DUSTY for the 
latter). Both spectra are individually rotationally and instrumental-Gaussian broadened, and then coadded 
in the ratio of the adopted spot covering fraction (for a covering fration /, the final spectrum is given by / 
times the spot spectrum + (1 — /) times the photospheric spectrum). The addition is performed prior to 
the normalization described above, to preserve the ratio of the spot to photospheric flux arising from their 
differences in both temperature and covering fraction. 

Finally, we note that the synthetic spectra were originally constructed at intervals of 100K and 0.5 dex 
in T c g and logg respectively. We have linearly interpolated between adjacent spectra (before normalization) 
to construct a finer final grid of models, with steps of 50K in T e s and 0.25 dex in \ogg. 



5. Results 

We begin with a general presentation of the fitting results and then quantify the uncertainties in the 
best fits. 

TiO-e: Fig. [T] shows the comparison between data and synthetic spectra in the TiO-e region, which has three 
bandheads at ~ AAA8435, 8445, 8455. We plot models at logg = 3.0, 3.5 and 4.0 (bracketing the empirically 
known value of 3.5), and T e g = 2400-2700K in steps of 100K (2650K is also shown to facilitate a comparison 
to K I, as we will discuss shortly). As expected, the model TiO is hardly sensitive to gravity over the 1 dex 
range plotted, but highly sensitive to temperature, with the bandheads at AA8445,8455 rapidly strengthening 
with decreasing T e g . We see that the 2500K model (in red) clearly fits the data very well, while cooler and 
hotter models (in blue) just as clearly do not. Given that even a 100K deviation from the best fit is evident 
to the eye, our precision in T e ff determination by eye is likely to be ~ ±50K (in agreement with M04). From 
the TiO-e fits, therefore, we would infer T cS w 2500K±50K. 

K I and TiO-e: Fig. [2] shows the comparison between the observed and model K I A7700, over the same 



5 'Pseudo-' because there is no true continuum in such cool objects, only an apparent continuum made up of millions of 
overlapping molecular lines; in the interests of conciseness, we drop the 'pseudo-' appelation forthwith. 
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range of log y and T c s plotted in Fig. [T]for TiO-e. We immediately see, as pointed out by |M04j, that the 
synthetic K I is sensitive to both gravity and temperature, becoming rapidly broader and deeper with either 
increasing gravity or decreasing T e g. At the T e g = 2500K inferred from TiO-e above, the best fit to K I is 
clearly obtained at logy=3.0. The model at logy — 3.5, for the same T e g, is obviously discre pant w ith the 



data; a conservative estimate of the error in this by-eye fit is ~±0.25 dex (in agreement with |M04[ ) . Thus, 
simultaneous model fits to TiO-e and K I imply T c g « 2500K and logy 3.0. This inferred gravity is lower 
than the empirical value of 3.52±0.03 by ~0.5 dex. 

Fig. [5] also shows that, if we impose the empirical log g of 3.5, a very good fit to K I is obtained at T e g » 
2650K. This is ^150K higher than derived from the TiO-e fits alone; as Fig.[T]shows, the synthetic TiO-e are 
incompatible with this temperature. We note that 2650K is consistent with T e g estimates for fiel d dwarfs of 



the same ~M6.5 spectral type as the primary (e.g.. iGolimowski et al.ll2004t ISlesnick et al.ll2004f l . while the 



2500K indicated by TiO-e is somewhat low in comparison. 

To better quantify the uncertainties in the fitting results of Figs. 1 and 2, we show the chi-square results 
for the data- model comparisons in Fig. 3. The TiO comparisons were carried out over the wavelength range 
[8420:8480] A (which includes all three bandheads; see Fig. 1), while the KI comparisons were over the 
range [7699: 7706] A (corresponding to the entire line, till the pseudo-continuum is reached on either side of 
line-center; see Fig. 2). These ranges correspond to 69 data points for KI, and 599 for TiO. 

The plot clearly shows that a degeneracy exists between T e g and log g in the KI line, as discussed above. 
There are two global minima over the range of temperatures and gravities examined, one at T e g sa 2650K 
and logy ~ 3.5, and another at T e g » 2850K and log g « 4.0. That is, a ^200K decrease in temperature 
can compensate for a 0.5 dex decrease in gravity (as also found by M04). Thus, at the empirical logy of 
2M0535, (logy = 3.5), the best-fit T e g from KI is 2650K, in agreement with our by-eye estimate above. The 
corresponding 1-er uncertainty is 30K. The TiO lines strongly indicate T e g = 2500±10K, again as found 
above. At this T e g, the K I line is marginally well fit at the 3-er contour level with logy = 3.0±0.05. Thus, 
the detailed chi-square comparisons are in excellent agreement with our fits-by-eye results for T e g and logy. 

To summarize: (1) fits to TiO-e yield T cS = 2500±10K; (2) adopting this T cB , fits to K I A7700 imply 
logy = 3.0±0.05, which is ^0.5 dex lower than the known gravity of the primary; and (3) adopting the 
known value of logy = 3.5 instead, fits to K I imply T c g = 2650±30K, which is ~150K higher than, and 
incompatible with, the value from TiO-e alone. 

We note here that, while the best-fit chi-square values are the same as obtained by eye, the formal 
uncertainties cited above for the chi-square analysis, obtained via interpolation over the model grid, are 
significantly smaller than the model grid spacing (50K in T c g and 0.25 dex in logy). We therefore adopt the 
grid spacing of 50K and 0.25 dex as a conservative estimate of our uncertainties for the rest of the paper 
(corresponding to the same errors assumed for the fits by eye). 

The discrepancies embodied in the above results may be due to lacks in the synthetic spectra, or an 
indication of real photospheric conditions. We discuss each in turn below. 



6. Discussion: Possible Interpretations of the Discrepancies in T e g and logy 
6.1. Model Opacity Uncertainties 

Reinersl |2005l hereafter R05) has compared the synthetic TiO spectra to observations of a sample of 
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early to mid-M field dwarfs, whose T e g and surface gravities are well-constrained via interferometric radius 
measurements. He shows that the model TiO-e bands systematically underestimate the temperatures of 
thes e objects. Assuming that uncertainties in the e-band model oscillator strengths — / e i(e) — are to blame, 
R05 estimates that / e i(e) 70% higher than adopted in the models would remove the discrepancy in the field 
dwarfs. Moreo ver, problems with / e i(e) should produce an analogous effect in young low-gravity brown 
dwarfs as well. Ir03 predicts that a 70% underestimation of / e i(e) would yield a 150-200K underestimation 
of T c ff in such young dwarfs, and an attendant \ogg (derived by imposing this T e g on the gravity-sensitive 
alkali lines, as in our analysis above) too low by ~0.3 dex. These results are in qualitative and quantitative 
agreement with the results presented above for 2M0535. 

The above prediction for young cool dwarfs is predicated on the assumption of problematic oscillator 
strengths. This is by no means proved, however; it may be that the fault lies in the adopted model equation 
of state instead (A. Reiners and P. Hauschildt, private comm., 2009), perhaps related to uncertainties in 
dust formation. In the latter case, it is not evident that the field dwarf discrepancies would necessarily be 
replicated in young objects. The steps required to resolve this issue are discussed further in Sj7] For now, we 
simply point out that the direction and magnitude of the / e i(e)-related T e g and log g discrepancies predicted 
by R05 for young brown dwarfs are consistent with those found in our analysis of the 2M0535 primary. As 
such, lacks in the synthetic spectra remain a viable explanation for our results. 

6.2. Metallicity Effects 

Alternatively, we have assumed a solar metallicity for the 2M0535 system (i.e., used synthetic spec- 
tra with log[M/H] = log[(M/H)/(M©/H0)] = 0.0); can non-solar abundances resolve the T e s and \ogg 
discrepancies we find? We do not think so, for the following reason. 

Metallicity variations affect the spectra as follows. Higher metallicity reduces the number of Hydrogen 
particles (which are the main source of collisional broadening) relative to metals; it also implies a decrease 
in pressure at a given optical depth (because of higher opacity) . Both effects tend to yield a narrower alkali 
line at higher metallicity, just as decreasing gravity does at solar metallicity (Schweitzer et al. 1996; Basri et 
al. 2000; Mohanty et al. 2004). This could lead to an underestimated log g from K I. Simultaneously, higher 
metallicity leads to an increase in the relative abundance of Titanium and Oxygen, and also causes a decrease 
in temperature at a given optical depth (again due to higher opacity) . Both effects lead to a strengthening of 
the TiO bands at higher metallicity, just as decreasing T e g does at solar metallicity (Leggett 1992; Mohanty 
et al. 2004). In summary, a higher metallicity mimics lower gravity and lower T c g. Thus, if we have 
underestimated the metallicity, we will also erroneously underestimate T e g and log g to compensate (i.e., to 
match the observed line profiles). 

These effects may be potentially invoked to explain our results in the following way. If metallicity in 
the 2M0535 system is higher than solar, then accounting for this will produce a T e g (from the TiO bands) 
that is somewhat higher than we currently find assuming solar abundances. Simultaneously, using the 
putative, higher-than-solar abundance would also lead to a higher log g (at any chosen T e g), from the KI 
line analysis, than we find at present. If the metallicity were sufficiently higher than solar, then these trends 
could potentially lead to an agreement in T e s between the TiO and KI regions at the correct (empirically 
known) gravity, resolving the discrepancies in our present results. 

However, it is the magnitude of the metallicity change required that is the stumbling block. On the 
one hand, Padgett (1996) has analyzed a number of nearby star-forming regions (Taurus- Auriga, Ophiuchus, 
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Chameleon and Orion, the latter being the region of which 2M0535 is a member) and found a solar abundance 
to within ±0.1 dex (i.e., ±25%) in all of them; within a given region, the variation is also at most 0.1 dex. 
On the other hand, Schweitzer ct al. (1996) have shown that, for a fixed T e g (~ 2700K, i.e., around the 
T e g regime of interest here), a 0.5 dex increase (decrease) in log[M/H] mimics a 0.5 dex decrease (increase) 
in logg in the synthetic alkali line profiles. Similarly, from our synthetic spectra, for a fixed log g we find 
that a 0.5 dex increase (decrease) in metallicity mimics a ^100K decrease (increase) in T e g in the synthetic 
e-band TiO @. 

Thus, the 0.1 dex maximum observed variation in metallicity would lead to non-significant changes in our 
inferred T c g and logg: ~20K and 0.1 dex respectively. These deviations are less than or comparable to the 
error bars on our derived values, and more importantly, completely insufficient to explain the discrepancies 
in temperature and gravity we find between the TiO and KI regions. We therefore posit that metallicity 
variations are not likely to explain our results, as an improbably large [M/H] > 0.5 dex for 2M0535 would 
be required. 



6.3. Cool Spots 



Finally, as discussed in m lCGB07l predict a significant presence of cool spots on the primary, to explain 



its T e ff reversal compared to the secondary. Assuming an admittedly extreme spot temperature of K, 
they require a spot covering fraction of ~50% (combined with severe magnetically-induced suppression of 
interior convection) to replicate the observations. More recently, a dopti ng more realistic spots 10% cooler 
than the bare photosphere in their photometric lightcurve analysis, find that a spot coverage of ~65% 



can reproduce the observed T e ff suppression of the primary (though, as mentioned in |JTJ these spots must 
be distributed axisymmetrically, to remain consistent with the relatively small photometric variations G09l 
observe; the latter require only <10% coverage by ncm-axisymmetric spots). 

In light of this, we investigate the effects of cool spots on our spectra. A priori, the following trends 
are expected. First, since the spots are by definition cooler than the unspotted photosphere, the TiO 
bandheads arising inside them will be deeper (relative to the continuum) than those from the surrounding 
photosphere; the resultant average TiO in a spatially unresolved spectrum will then imply a temperature 
intermediate between that of the spotted and unspotted surfaces. This trend will not however be monotonic 
with decreasing spot temperature, since the continuum flux from spots also falls with the spot temperature: 
the spot TiO contributes less to the average TiO as the spot gets cooler. Thus, for a fixed unspotted 
photospheric T e g and spot coverage, the average TiO in the combined spotted±unspotted spectrum first 
deepens (relative to the TiO from an unspotted surface) with decreasing spot temperature, and then reverses 
as the spots become still cooler, to become shallower again (i.e., approaches again the TiO from an unspotted 
surface). The spot temperature at which this reversal occurs, and thus the maximum depth of the TiO 
bandheads in the presence of spots, is determined by both the unspotted photospheric T e ff and the spot 
coverage assumed. As an extreme example, K spots contribute no flux at all, and hence, for any coverage 
<100%, will have no effect on the shape of the TiO, which arises in this case only from the unspotted surface 
(even though the absolute flux in the continuum and bandheads will be lower than in the absence of spots, 
since the unspotted surface now covers only a fraction of the total stellar surface). The trend in TiO with 



6 Only low-resolution, logg = 5.0 synthetic spectra are currently available for non-solar metallicities in the COND and DUSTY 
models. However, the T e ff versus log[M/H] changes are clear even at low-resolution, and the insensitivity of the TiO e-band to 
log g indicates that these results should be reasonably applicable at the lower gravity of our sources as well. 
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changing spot temperature and coverage is illustrated in the top panel of Fig. 4. 

Second, inside a spot, magnetic fields provide partial support against the external photospheric gas 
pressure; co nsequently, the gas pressu re alone within a spot is lower than in the surrounding unspotted 



surface (e.g., lAmado et al.lll999l l2000h ■ This exactly mimics the reduction of gas pressure caused by lower 
surface gravity, and causes the gravity- (more accurately, gas-pressure-) sensitive alkali absorption lines (e.g. 
K I) in a cool spot to be narrower than outside, thereby implying a lower effective gravity within the spot. 
The averaged alkali lines in the combined spotted-l- unspotted spectrum will then imply a gravity intermediate 
between the true photospheric value and the effective one in the spot. Conversely, spots are also cooler than 
the external photosphere; this tends to make the alkali lines, which are temperature-sensitive as well, broader 
within a spot, mitigating the narrowing caused by the magnetic pressure effects and reducing the apparent 
gravity offset. Finally, both effects are limited by the decreasing flux from a spot with lower temperature, 
exactly as discussed above for TiO; for a given unspotted photospheric T e ff and spot coverage, the shape 
of an alkali absorption line is negligibly changed, relative to that from an unspotted surface, once the spot 
temperature falls below a certain threshold. The trend in KI with changing spot temperature and coverage 
is illustrated in the bottom panel of Fig. 4. 

These trends have the following consequences. As we have shown, the maximum depth of the TiO 
bandheads in the presence of spots is determined by both the unspotted photospheric T c g and the spot 
coverage adopted. Thus, for an observed TiO band depth, specifying the unspotted photospheric T e g sets a 
lower limit to the spot areal coverage: if the coverage is below this value, the average TiO band depth reverses 
before it can match the observed depth regardless of how cold the spots become. For a coverage higher than 
this threshold and fixed unspotted photospheric T e s, there is a degeneracy between spot temperature and 
coverage: a larger coverage allows cooler spots. At the same time, if we decrease the spot temperature we 
must also decrease the spot's effective gravity in order to match the observed alkali absorption lines, given 
the competing effects of gravity and temperature in these lines. 

These trends imply that, with a priori knowledge of the unspotted T ff , and the true stellar gravity as 
well as the effective gravity within a spot, one can solve for the spot temperature as well as covering fraction 
via simultaneous analysis of the TiO and alkali absorptions. We however only know the stellar gravity, and 
have no advance knowledge of either the unspotted photospheric T c ff or the spots' effective gravity. As such, 
we can make no claims to a unique solution, or to a full search of the available parameter space. Instead, 
our goal is a plausible solution, based on constraints set by known properties of sunspots and starspots in 
general as well as by the 2M0535 lightcurve analysis so far. In pa rticular, we assume that the spots are 



cooler than the surrounding photosphere by at most ^10-25% (e.g.. lG09i ; iLinskv et al.ll2002l ). and that the 
differential between their effective gravity and the higher photospheric value (logg=3.5) is < 0.5 dex (e.g., 
Amado et al1ll999l . lioooh . The unspotted photospheric T c ff is also kept a free parameter, but with a lower 



bound of 2500K set by the T e g inferred from the TiO fits in §5 (the reason for this lower bound is that, if the 
unspotted photosphere T e g were < 2500K, then combining it with spots that are by definition even cooler 
could never produce TiO bands that appear to be at 2500K when compared to unspotted models). 

Within these constraints, we construct star+spot models by assigning synthetic spectra to the unspotted 
and spotted surfaces, as detailed in g] We find a viable solution for our spectroscopic data with the 
following parameters: an unspotted stellar surface with T e g = 2700K and (empirically determined) log g = 
3.5, combined with 70% axisymmetric areal coverage by spots with a temperature of 2300K and effective 
logy = 3.0 (i.e., 15% cooler, and 0.5 dex lower apparent gravity, than the unspotted surface). 



Fig. 5 demonstrates the viability of this spotted model. In the left panels, we compare the data (in 
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black) to the best-fit unspotted models from Figs. 1 and 2 (TiO: best-fit T G ff of 2500K [in red\, as well as 
2550K [in grey) to show the deviation caused by our adopted 50K uncertainty; KI: best-fit log g = 3.0 at 
the 2500K implied by TiO [in blue], as well as best-fit T e s = 2650K at the empirically determined log g of 
3.5 [in red\). In the right panels, we compare the best- fit spot model described above (in green) to both the 
data and the best-fit unspotted models from the left panels. We see that the spotted model very closely 
reproduces the data as well as the unspotted model fits (in particular, the spotted model is identical to the 
2550K unspotted model in TiO, i.e., within 50K - our adopted error - of the best-fit 2500K unspotted TiO 
model; it is also nearly indistinguishable from the 2650K KI model at the empirical gravity of log g = 3.5). 

Thus, the discrepancies in T c r and logg implied by the unspotted model-fits to TiO-e and K I are 
resolved by this single star+spot model. Consequently, cool spots are a viable explanation for our data. 



7. Summary and Conclusions 

We have shown that the TiO-e and K I absorption features in the high-resolution optical spectrum 
of the higher-mass primary in 2M0535 (2M0535A) are consistent with a T e g = 2700K, (empirical) logg = 
3.5 photosphere combined with cool spots with a temperature of 2300K and effective gravity of log q = 
3.0 covering 70% of the brown dwarf's surface. This is in agreement with the scenario outlined bv lCGB07l . 
wherein the temperature reversal in the primary relative to the secondary is caused by magnetic fields, which 
induce both a reduction in the convective efficiency and hig h cool spot coverage. While the ext reme spot 



covering fraction we find is similar to that inferred bv lCGB07l (50% using unrealistic K spots) and lG09l (65% 



with more realistic spots 10% cooler than the photosphere), this very high fraction is nevertheless troubling: 
in effect, it makes 2M0535A appear to be a "very cool" (2300K) star covered by hot spots, rather than the 
reverse. On the other hand, we note that Stauffer et al. (2003) argue that the anomalous colours of Pleiades 
K and M dwarfs result from axisymmnetrically distributed cool spots with a very large areal coverage, > 
50%; Gullbring et al. (1998) have argued for similarly large axisymmetric spots, with covering fractions of 
~ 50-70%, to account for the anomalous colors of even younger weak-lined T Tauri stars (WTTS). There is 
also evidence for the presence of such axisymmetric large spots from Doppler imaging of WTTS, e.g., large 
polar spots in V410 Tau (Hatzes 1995) and HDE 283572 (Joncour et al. 1994). Thus, such spots may be 
usual during the early evolution of these stars, when they are rapidly rotating and highly active, and the 
phenomenon may extend into the substellar regime as well. 

The other option, which we show our results are also consistent with, is that the synthetic spectra are 
in error, leading to a discrepancy between the T e g and log g deri ved from simultaneous fits to TiO-e and K I 
(which then leads us to postulate a surfeit of cool spots). R05l postulated such synthetic spectr 



trum errors, 



also seen in analyses of field dwarfs, to arise from problems with the model TiO-e oscillator strengths. As we 
were submitting this paper, it came to our attention that the model errors may lie in the adopted equation of 
state instead, and that newer models rectifying this are being prepared (A. Reiners & P. Hauschildt, private 
comm., 2009). Whether these models can resolve the discrepant values obtained from TiO-e and K I in the 
case of the very young 2M0535A, without having to resort to copious cool spots, remains to be seen. 

Regardless of whether the new models fare better or not, setting rigorous constraints on both the models 
and the physical conditions on 2M0535A — i.e., determining whether the model fits (and thus the implied T c s 
and effective log g for the spots) are truly valid and/or if very large cool spots exist on the primary — now 
requires independently carrying out exactly the same analysis for the lower-mass secondary (2M0535B). 



CGB07J's theory predicts that the secondary should be much less spotted than 2M0535A: thus, if the same 
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discrepancies between TiO-e and K I appear in the secondary as well, then errors in the models will be 
clearly implied; if not, then the suggestion that the primary has an extremely large spot covering fraction 
will be bolstered. We are currently undertaking observations of 2M0535B to carry out this test. 

Finally, if it turns out from 2M0535B's analysis that the current models are in error (i.e., if the same 
discrepancies between TiO and K I appear in 2M0525B as well even though it is expected to be relatively 
unspotted) , and thus there is no observational rationale for suggesting very large cool spots on the primary 
component 2M0535A, then the cause of the temperature reversal between the two binary components once 
again becomes an unresolve d issue. In that case, it may be that the theory proposed very recently by 
MacDonald fe Mullanl (|2009T ) — wherein magnetic fields are again to blame, but by inhibiting the onset of 
convection (though not complete ly) throughout the star, instead of just in the uppermost super-adiabatic 
layers as in the theory of ICGB07I — may be the correct one. Again, we stress that this can only be tested via 
comparison with an analogous analysis of 2M0535B. Finally, complementary high-resolution spectroscopic 
obse rvations in the near-infrared (NIR) can be very useful in determining the properties of cool brown dwarfs 
(e.g. Zapatero Osorio et al. 20061 : McLean et al. 2007 : Del Burgo et al. 200!):. as brown dwarfs are brightest 
in the NIR. This would be particularly useful for a comparative analysis for the 2M0535 system, where there 
is a more favorable contrast ratio between the secondary and the primary in the NIR; it would also be helpful 
for assessing the importance of cool spots, whose influence is less marked in the NIR than in the optical. 

This binary, the first eclipsing system in the substellar domain, has already proved to be a rich testing 
ground for theories of brown dwarf formation and evolution, and it promises to continue being a Rosetta 
Stone in this regard. 
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acknowledges support from NSF grant AST-0607773. This material is based upon work supported by AURA 
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Fig. 1— Observed TiO-e region in 2M0535A (black) compared to Dusty (< 2500K) and Cond (> 2500K) 
models. Best-fit model (Cond 2500K) shown in red; all others, which clearly diverge from the data by eye, 
shown in blue. Note that the model fits are very insensitive to gravity over the 1 dex range plotted. See §5. 
This figure is shown in color in the electronic version only. 
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Fig. 2.— Observed red lobe of K I in 2M0535A (black) compared to Dusty (< 2500K) and Cond (> 2500K) 
models. Best-fit models (Cond [2500K, log g = 3.0] and [2650K, log g = 3.5]) shown in red; worse but still 
admissible fits by eye shown in purple; and all others, which clearly diverge from the data, shown in blue. 
Note that the K I absorption is sensitive to both T e g and gravity: a 150K increase in T e g compensates for a 
0.5dex rise in log g. At T e g = 2500K, corresponding to the best-fit to TiO-e (Fig. 1), the K I implies log g 
= 3.0, while at the empirically determined log g = 3.5, it implies T e g = 2650K. See §5. This figure is shown 
in color in the electronic version only. 
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Fig. 3. — Determination of goodness-of-fit and formal fit parameter uncertainties. Contours of constant 
X 2 — Xmin = 2.3,6.2,11.8, representing 1,2,3-ct joint confidence intervals in the T c g~\ogg parameter plane. 
Top: Joint confidence intervals for fitting of K I. The absolute minimum x 2 best fit is for T c tf = 2850 K and 
\ogg = 4.0, however a second equally good fit within 1-a confidence occurs at T e g = 2650 K and log g = 3.5. 
Bottom: Joint confidence intervals for fitting of TiO. The contours demonstrate that for TiO the best fitting 
model spectra are relatively insensitive to log g but highly sensitive to T e s; a best- fit T e g = 2500 ± 50 is 
strongly preferred at high confidence. 
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Fig. 4. — Trends in TiO (top panel) and KI (bottom panel) for spotted photospheres with varying spot 
temperature and areal coverage. For all models in both panels, the unspotted photosphere has T e g — 2700K 
and log g = 3.5, while the spots have an effective gravity of 3.0. In both panels, the black solid line shows the 
unspotted model; the red solid line shows the model with [Tspot, coverage] = [2300K, 50%]; the blue solid 
line shows the model with [Tspot, coverage] = [2100K, 50%]; and the blue dashed line shows the model with 
[Tspot, coverage] = [2100K, 85%]. We see that the [2300K, 50%] spot model is deeper in TiO and broader 
in Ki than the unspotted case, but the 2100K spot model with the same areal coverage of 50% has reversed 
in strength, and is closer to the unspotted case than the 2300K model. Increasing the areal coverage of the 
2100K case to 85%, however, makes it significantly deeper in TiO and broader in KI than both the unspotted 
and 2300K spotted models. These changes illustrate the spot-related trends discussed in §6.3. This figure is 
shown in color in the electronic version only. 
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Fig. 5. — Data for 2M0535A (black) compared to unspotted models (left panels) and a star+spot model 
(right panels). The TiO comparisons are shown at the top and the K I comparisons at the bottom. The 
unspotted models plotted are the same best-fit ones plotted in Figs. 1 and 2 (one best-fit model for TiO (in 
red), two for K I (in red and blue)). Additionally, we plot a 2550K model (in grey) for TiO, to illustrate our 
adopted error bar of 50K in the T e ff derived from TiO. The spotted model (in green in the right panels) is a 
[2700K, log g = 3.5] photosphere with [2300K, log g = 3.0] spots covering 70% of the surface. The fits to the 
spotted model are nearly indistinguishable from the fits to the unspotted models (specifically, the spotted 
model is an excellent match to the 2550K unspotted model in TiO, i.e., within our adopted T c ff error bar 
for TiO, and to the [2650K, logg = 3.5] unspotted model in KI), implying that the spotted model is as good 
a description of the data as the unspotted ones. See §6.3. This figure is shown in color in the electronic 
version only. 



